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LONG-TERM  GOALS 

Our  main  goal  is  to  assess  the  hydrodynamics  of  turbidity  currents  and  mudflows  and  their  ability  to 
produce  morphological  features  such  as  ripples,  dunes,  antidunes  and  gullies  along  their  paths.  Of 
particular  interest  is  to  elucidate  the  role  played  by  turbidity  currents  on  the  inception  of  submarine 
channels  and  canyons.  To  this  end,  understanding  the  mechanics  of  sediment  erosion,  including 
bedrock,  by  density  currents  is  another  long-term  goal  of  this  research. 

OBJECTIVES 

Our  main  objective  is  to  advance  our  understanding  of  the  seabed  morphology  associated  with  the 
passage  of  turbidity  currents  and  mudflows.  Different  morphological  features  encountered  in  the 
marine  environment  are  related  to  the  action  of  bottom  currents  and  turbidity  currents.  Large  fields  of 
long-wavelength,  upslope-migrating  bedforms  observed  on  submarine  channel  levees  are  attributed  to 
the  action  of  unconfined  turbidity  currents  that  overflow  channels  and  canyons  and  spread  sediments 
across  levees.  Although  it  has  been  largely  argued  that  such  sediment  waves,  with  ripples 
superimposed  on  them,  are  the  result  of  downslope-flowing  density  underflows,  the  origin  of  such 
features  still  remains  under  debate. 

APPROACH 

Our  approach  consists  of  a  combination  of  theoretical  analysis,  numerical  modeling,  and  laboratory 
experiments.  Stability  analysis  with  multiple  scales  was  used  to  investigate  whether  supercritical 
conservative  density  currents  are  capable  of  developing  long-wavelength,  upstream-migrating 
bedforms.  Laboratory  experiments  were  performed  in  the  MARGINS  TANK  to  observe  if 
conservative  density  currents  and  depositional  turbidity  currents  are  capable  of  developing  bedforms  at 
different  scales  such  as  dunes  and  ripples,  as  well  as  long-wavelength,  antidune-like,  sediment  waves. 
Two-  and  three-dimensional  direct  numerical  simulations  of  a  discontinuous  density  current 
propagating  down  a  slope  were  conducted  to  assess  the  role  of  mixing  on  the  flow  dynamics. 

WORK  COMPLETED 

A  theoretical  stability  analysis  has  been  completed  using  conservative  density  currents  capable  of 
transporting  sediment  as  bedload  as  well  as  erosional  turbidity  currents.  Comparison  of  maximum 
growth  rate  curves  for  the  above  mentioned  cases,  as  function  of  Richardson  number,  show  that 
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turbidity  currents  tend  to  fonn  longer  sediment  waves  than  those  associated  with  conservative  density 
currents  with  bedload  transport  only.  The  analysis  has  been  extended  to  account  for  the  role  played  by 
the  flow  turbulence  in  the  development  of  bedforms,  by  adding  a  conservation  equation  for  turbulent 
kinetic  energy. 

More  than  thirty  laboratory  experiments  were  completed  in  the  MARGINS  TANK  for:  a)  wide 
turbidity  currents  emanating  from  a  line  source  and  producing  incipient  gullies,  and  b)  channelized 
turbidity  currents  generating  a  range  of  bedforms  from  ripples  to  dunes  and  antidunes. 

Several  two-  and  three-dimensional  numerical  simulations  of  discontinuous  density  currents  were 
conducted.  The  numerical  results  were  compared  against  unpublished  laboratory  observations 
previously  made  by  the  PI  at  the  University  of  Minnesota.  The  numerical  simulations  were  carried  out 
using  a  stabilized,  equal-order  finite  element  method,  in  which  a  discontinuity-capturing  technique  was 
embedded  for  the  salt  transport.  This  technique  captures  sharp  fronts  with  little  spurious  oscillations, 
allowing  local  features  such  as  shear  instabilities  at  density  interfaces  to  be  well  reproduced. 

RESULTS 

Results  obtained  from  the  bed  stability  analysis  indicate  that  both  conservative  and  non-conservative 
supercritical  dense  underflows  are  able  to  develop  long-sediment  waves  that  migrate  upslope, 
resembling  antidunes  (Fedele  and  Garcia,  2005a  and  2005b).  Figure  1  shows  the  components  of  the 
complex  celerity  of  bed  perturbations,  as  a  function  of  wave-number  (K),  for  a  density  current.  In  this 
figure,  growth  rate  (C/  >0)  indicates  that  long  sediment  waves  are  likely  to  develop,  whereas  the 
celerity  (Cr  <  0)  indicates  an  upstream  migration  for  these  bed  instabilities.  One  of  the  important 
results  obtained  from  the  stability  analysis  is  the  spatial  dependence  of  the  (complex)  dispersion 
relationship,  which  gives  the  growth  rate  of  bed  perturbations.  Figure  2  shows  growth  rates  of  bed 
perturbations  computed  using  the  imaginary  component  of  the  dispersion  relationship,  for  two  different 
locations  along  the  flow  path.  It  is  observed  that  antidunes  formed  by  supercritical  conservative 
density  currents,  with  bedload  transport  only,  tend  to  elongate,  as  the  maximum  growth  rate  moves 
towards  smaller  wavenumbers  in  the  downslope  direction.  Also,  the  value  of  the  positive  peak  of  the 
growth  rate  tends  to  decrease  in  that  direction  as  well,  indicating  that  wave  amplitude  tends  to  decrease 
as  the  current  flows  downslope. 

Laboratory  experiments  showed  that  conservative  density  currents  and  depositional  turbidity  currents 
are  capable  of  developing  bedfonns  at  different  scales  such  as  dunes  and  ripples,  as  well  as  long- 
wavelength  antidune-like  sediment  waves  (Fedele,  2004).  The  laboratory  observations  in  the 
MARGINS  tank  indicated  that  the  front  of  turbidity  currents  plays  an  important  on  the  bottom 
morphology.  In  particular,  in  order  to  be  able  to  model  the  inception  and  development  of  gullies,  it  is 
necessary  to  have  a  3-D  hydrodynamic  model  of  the  flow. 

NUMERICAL  RESULTS 

A  rigorous  model  for  two-way  coupled  multiphase  flow  has  been  developed.  The  model  has  been 
derived  formally  from  the  complete  two-phase  flow  models  presented  in  Zhang  and  Prosperetti  (1997) 
and  Balachandar  and  Ferry  (2004).  In  contrast  to  the  models  presented  by  these  authors,  the  model  that 
we  have  developed  is  based  on  the  volume-averaged  velocity  instead  of  the  continuous  phase  velocity. 
In  this  way,  by  imposing  the  divergence-free  condition  to  the  velocity  field,  we  satisfy  exactly  mass 
conservation  since  for  incompressible  phases  mass  conservation  reduces  to  volume  conservation.  This 
feature  of  the  model  facilitates  its  implementation  in  standard  incompressible  Navier-Stokes  solvers. 
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The  derived  momentum  conservation  equations  present  three  terms  in  addition  to  the  standard 
equations  that  account  for  the  two-way  coupling  between  phases.  The  model  is  completed  by  a 
transport  equation  for  the  void  fraction,  which  is  transported  by  the  disperse  phase  velocity.  For  the 
disperse  phase  velocity  we  have  derived  a  new  version  of  the  equilibrium  Eulerian  method  based  on 
the  volume-averaged  velocity.  The  derivation  of  the  model  was  done  by  expanding  all  the  terms  using 
the  small  parameters  x  and  a  (particle  response  time  and  volume  fraction,  respectively)  and  retaining 
tenns  0(1),  0(x)  and  0(a),  which  makes  the  model  exact  to  0(xa+x  +a").  Finally,  it  is  worth 
mentioning  that  the  model  is  the  most  sophisticated  version  of  the  equilibrium  Eulerian  method:  its 
innovations  include  exact  multiphase  volume  conservation,  finite  settling  velocity  effects,  and  the 
incorporation  of  volume  fraction  gradient  and  Faxen  terms.  The  derived  model  has  been  implemented 
in  a  spectral  DNS  channel  flow  code,  which  allows  for  the  detailed  resolution  of  fine-scale  features 
without  resorting  to  subscale  modeling. 

We  have  perfonned  3D  simulations  for  planar  and  cylindrical  configurations  (Cantero  et  al.  2004a, 
Cantero  et  al.  2005c)  with  Gr=1.5  x  106  and  Gr=1.5  x  107.  The  results  are  in  figures  3  and  4.  These 
solutions  not  only  shows  the  Kevin-Helmholtz  billows,  but  also  the  existence  of  a  fast  instability  at  the 
front  of  the  current  near  the  vertical  boundaries  that  develops  into  a  complex  lobe  and  clef  pattern.  We 
can  also  observe  the  existence  of  an  instability  that  develops  at  the  interface  of  the  current.  Figure  5 
shows  the  time  evolution  of  the  lobe  and  cleft  instability.  In  this  figure  the  front  is  visualized  by 
contours  of  density.  Figure  6  shows  the  radius  of  curvature  of  the  front.  Dark  regions  correspond  to 
negative  curvatures  that  are  characteristic  of  clefts.  Light  color  regions  correspond  to  positive 
curvatures  that  are  characteristic  of  lobes.  Studies  generalizing  these  results  to  particulate  turbidity 
currents  are  under  way  (Cantero  et  al.  2005b). 

We  have  also  performed  studies  on  the  particles  influence  on  the  front  velocity  of  gravity  currents 
(Cantero  et  al.  2004b,  Cantero  et  al.  2005a).  Figure  7  shows  the  results  for  the  limit  when  particles  are 
so  small  that  they  act  as  a  scalar  field  (x  =0  and  w=0).  The  flow  is  visualized  by  a  contour  of 
dimensionless  particle  concentration  =0.1.  Soon  after  the  release  an  intrusion  front  forms  with  a  lifted 
nose  due  to  the  no-slip  boundary  condition.  As  the  current  advances,  Kelvin-Helmholtz  vortices  are 
shed  from  the  front  which  produce  a  net  drag  that  balances  the  initial  acceleration  of  the  front.  As  a 
consequence,  after  the  initial  set-up  of  the  Kelvin-Helmholtz  vortices,  the  front  moves  at  constant 
speed  until  the  dilution  in  the  current  becomes  important.  Then,  the  current  slows  down  and  eventually 
dissipates. 

Figures  8  and  9  show  the  results  for  currents  of  inertial  particles  with  negligible  settling  (x=0.05  and 
w=0  and  x  =0.1  and  w=0).  Two  contours  are  shown  in  these  figures:  the  solid  line  contour  that 
corresponds  to  particle  concentration^.  1,  and  the  dash  line  contour  that  corresponds  to  the  initial 
particle  concentration  equal  to  1.  Two  main  differences  are  observed  compared  to  the  no-inertia- 
particles  case.  The  first  one  is  the  migration  of  particles  away  from  the  core  of  the  Kelvin-Helmholtz 
vortices,  and  the  second  one  is  the  accumulation  of  particles  in  the  front  of  the  current,  producing 
regions  of  particle  concentration  greater  than  the  initial  value.  These  two  effects  can  be  explained  by 
noticing  that  the  divergence  of  the  particles  velocity  field  is 

v  •  u(<  =  f(i  -  (3)  (||n||2  -  ns  |2) 

where  S  and  Q  are  the  symmetric  and  skew-symmetric  parts  of  the  local  fluid  velocity  gradient  tensor 
and  p  is  the  density  ratio  of  the  phases.  Note  from  this  equation  that  for  particles  substantially  heavier 
than  the  continuous  phase  (P  =  0),  the  divergence  of  the  particle  velocity  field  has  the  sign 
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of  (| | W| |  -  ||S||  ),  which  means  that  heavy  particles  migrate  from  regions  of  vorticity  and  accumulate  in 
regions  of  high  strain  rate. 

The  preferential  particles  accumulation  described  above  has  a  very  important  consequence  in  the 
current  front  velocity.  Figure  10  shows  the  front  velocity  for  the  three  different  cases  mentioned  above. 
Observe  that  the  front  velocity  in  the  phase  of  constant  velocity  increases  for  larger  values  of  t.  The 
current  with  t=0.  I  presents  a  front  velocity  5%  larger  that  the  current  with  t=0. 

IMPACT/APPLICATION 

Theoretical  analyses  supported  by  laboratory  observations  have  clearly  shown  the  capability  of 
turbidity  currents  for  producing  longitudinal  small-scale  bedforms  such  as  ripples  as  well  as  long- 
wavelength  antidunes.  Our  understanding  of  the  mechanics  of  bedforms  in  continental  margins  until 
now  has  been  rather  limited  and  full  of  speculation.  We  hope  that  this  work  will  facilitate  both  the 
interpretation  of  the  geologic  record  as  well  as  the  design  and  placement  of  submarine  structures  on 
stable  sediment  deposits.  Our  work  suggests  that  turbidity  currents  could  also  be  a  mechanism 
important  for  mine  burial  and  should  therefore  receive  more  attention,  particularly  at  river  mouths 
where  high-sediment  discharges  can  generate  such  underflows. 

RELATED  PROJECTS 

We  are  collaborating  with  Pat  Wiberg,  University  of  Virginia,  and  Lincoln  Pratson,  Duke  University, 
to  develop  and  incorporate  a  field-tested,  mechanistic  model  for  sediment  erosion  into  our  density 
current  model  as  part  of  the  EuroStratafonn  Program. 
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FIGURE  1.  Stability  of  bed  perturbations  for  a  conservative  density  current  with  bedload  transport. 
Figure  shows  imaginary  and  real  components  of  the  celerity  of  bed  perturbations  as  functions  of 

bedform  wave  number  (x*=0). 


FIGURE  2.  Down  slope  variation  of  characteristic  wavelength  and  growth  rate. 
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FIGURE  3.  Solution  for  the  3D  planar  saline  density  current,  Gr=1.5  x  106.  From  top  to  bottom  t=0, 
t=5,  t=10  and  t=15.  The  surface  represents  a  void  fraction  isosurface. 
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FIGURE  4.  Solution  for  the  3D  cylindric  saline  density  current,  Gr=1.5  x  106.  From  top  to  bottom 
t=0,  t=4,  t=8  and  t=12.  The  surface  represents  a  void  fraction  isosurface. 
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FIGURE  6.  Time  evolution  of  curvature  of  lobe  and  cleft  instabilities.  Solution  for  the  3D 
cylindrical  saline  density  current,  Gr=1.5  x  106. 
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FIGURE  9.  Contours  of  particle  concentration  solid  line=  0.1  and  dashed  line=1.0.  Solution  for 

r=0. 1 ,  w=0  and  Gr=1.5  x  106. 


